Astronomy & Astrophysics manuscript no. skew'fl 


2nd February 2008 


(DOI: will be inserted by hand later) 





The cosmic shear three-point functions 



K. Benabed 1 and R. Scoccimarro 2 



in 
o 
o 



> 

00 
CN 

in 
o 
in 
o 

Oh 
i 

o 



1 Institut d'Astrophysique de Paris, 98bis bd Arago, 75014 Paris, France 

2 Center for Cosmology and Particule Physics, Department of Physics, New York University, New York, NY 
10003, USA. 

the date of receipt and acceptance should be inserted later 

Abstract We investigate the three-point functions of the weak lensing cosmic shear, using both analytic methods 
and numerical results from N-body simulations. The analytic model, an isolated dark matter halo with a power- 
law profile chosen to fit the effective index at the scale probed, can be used to understand the basic properties 
of the eight three-point functions observed in simulations. We use this model to construct a single three-point 
function estimator that "optimally" combines the eight three-point functions. This new estimator is an alternative 
to M a p statistics and provides up to a factor of two improvement in signal to noise compared to previously used 
combinations of cosmic shear three-point functions. 
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1. Introduction 

The quality of recent and upcoming galaxy weak lens- 
ing surveys is rapidly improving and allows increas- 
ingly high signal to noise determination of second and 
third order cosmic shear correlation functions, which con- 
tain very interesting cosmological information. For ex- 
ample, the two poi nt function constrains a combina- 
tion of Q, m and as l)Bernardeau et al1ll997|) with some 
sensiti vity to the shape of the primordial power spec- 
trum l|Schneider et alJl2002|) and the equation of state 
of the dark energy component ijBenabed fc Van Waerbekel 
l2004lk The three-point function, which we investigate in 
this paper, is particularly important for breaking degen- 
eracies on two-point statistics, giving a strong constraint 
on Q, m l|Bernardeau et alJll997j) with a small dependency 
on th e dark energy component l|Benabed fc Bernardeaul 
l200l|) . Studies of three-point statistics have recently re- 
ceived a great deal of attention from the theoretical 
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2003; Schneider & Lombardi 2003: 
120021 Zald arriaea fc Scoccimarrd 



20031 ISchneider et al.ll2005|) . inspired in part by the de- 
tection of an averaged c osmic shear three-point function 
( Bernardeau et al.ll200^l . 

Weak gravitational lensing detection is based 
on studying the alignment of background galaxies 
due to the lensing effect o f the intervening mass 
l|Bartelmann fc Schneider! l200l|) . which provides a mea- 
sure of the average shear 7. From this estimator, one 
could in principle obtain the projected mass n, however 
this requires inverting a non-local relation (see Eq. 



below) , which is very sensitive to the boundary conditions 
and thus difficult in galaxy surveys with complicated 
geometries and masks due to bright stars, etc. 

One possible alternative is to build statistics of the 
projected mass k without reconstruction of k itself, by 
finding a local operation on the shear that will give some 
filtered version of the projected mass. A well-known pro- 
cedure of this type, t he aperture mass statistic or M ap 
[ Schneider et al.ll2002|) . corresponds to averaging the two- 
point function of k with a compensated (zero integrated 
volume) filter. In terms of the shear, it corresponds to con- 
volving the shear two-point functions with a compact sup- 
port filter. This statistic is not optimal in terms of signal 
to noise, and although it can be extended to the three - 
point correlation function l|Schneider fc Lombardill2003|b 
it leads to a determination o f the projected skewness with 
a rat her poor signal to noise ijjarvis et alJl2004l IPen et aT] 
l2003h . This can be solved partly by suitably designing the 
window involved in the M ap statistic, but at th e price of 
more complicated equations (jjarvis et al.ll2004)l . 

For these reasons one would like to use statistics of the 
shear field to constrain cosmological models. The difficulty 
with this approach is that the shear is a spin-2 field, and 
thus the geometrical properties of its correlation functions 
are much more complicated than for the scalar field k, 
particularly when one goes beyond two-point statistics. 
For a general triangle configuration, there are eight non- 
vanishing three-point functions of the cosmi c shear field 
l|Scrmeider et a,l.ll2002HTa,kada fc .Talnll2003ah. 

In this work we study the main features of the cosmic 
shear three-point functions using numerical simulations 
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(see ijTakada fe Jainll2003aj) for previous work), and com- 
pare these r esults to a simple analytic mode l (following 
the ideas in i Zaldarriaea fc Scoccimarrol2003[) ') to extract 



the main features for the purpose of defining an optimal 
linear combination of the eight three-point functions into 
a single object, which is e asier to deal with. This g eneral- 
izes the previous study in l|Bernardeau et al . 2003), where 
a particular linear combination was selected based on the 
pattern of the three-point functions under some condi- 
tions, tested empirically with numerical simulations. Here 
we pay particular attention to avoiding cancellations that 
can occur when doing such linear combinations, with the 
purpose of maximizing signal to noise. 

This paper is organized as follows. In section [2] wc 
discuss the basics of weak lensing cosmic shear, and 
present the geo metrical prop erties of the shear two- 
point functions ijKaiserl ^992) an d then extend the re- 
sults to the three-point function s l)Schneider et al.l [2002: 
IZaldarriaga fc Scoccimarroll2003h . In section El we present 
a simple analytic model based on a power-law profile, 
which allows us to extract the basic properties of the shear 
three-point functions. In section 0] we present the results 
of measurements in numerical simulations and compare 
to the analytic model. In particular, we show that in the 
range of scales between 1 and 3 arc minutes they agree- 
ment with an isothermal sphere is very good. In sectional 
we use these results to construct an estimator that com- 
bines the eight three-point functions into a single object, 
and compare the signal to noise of this new statistic to 
the estimator previously used in t he literature to detect a 
cosm ic shear three-point function i Bernardeau et alT2 002. 
W 



2. Weak lensing and three point function 

2.1. The Convergence and Shear fields 

To first order in the density perturbations, the weak lens- 
ing effect is determined completely by its convergence 
field k which is a simple project ion of the matter density 
contrast along the line of sigh t ijBartelmann fc Schneider! 
l200lUBernardeau et al]ll997jl 

k{6) = j ^-w{z)8{d,z), (1) 

where the weight w is the lensing efficiency function which 
reads 

, .. 3„ D(z)D(z, z») , 
^ )= 2°° D(«.) ^ < 2 > 
in the case of a infinitely thin source plane at redshift 
z s . Here T)(z\,Z2) is the angular distances between red- 
shifts z\ and z-i and tt m the present dark matter density 
in units of the critical density. We will focus our anal- 
ysis on the three-point function of the lensing shear 7. 
Ignoring higher-order corrections perturbation theory, the 
shear field is simply related to the convergence field k by 
the non-local equations 

71 = A" 1 (df ~ dl) «, 

72 = A-^di&K, w 



where A -1 denotes the inverse Laplacian operator. With 
these equations, the information contained in the weak 
lensing effect is simple to understand. The convergence 
field is simply a projection along the line of sight of the 
matter density contrast. Thus its second and third mo- 
ments are related, up to projection effects, to the two- and 
three-point functions of the density contrast. The shear 
and convergence field, being equivalent through Eq. Q, 
contain the same information. 

The shear 7 transforms as a spin-2 object. Equation 
© is the analogous to that defining the E and B com- 
ponents of the polarization in terms of the Q and U 
Stockes parameters. The convergence field plays here the 
role of the E polarization, and there is no B compo- 
nent to the weak lensing effect to this order. Such con- 
tribution can be produced by deviations of the Born ap- 
proximation, and have been shown to be more than two 
orders of magnitude sma l ler than the do minant scalar 
mode (|Coorav fc HiJl2002l Ijain et alJl200(A . In addition, 
in observations B modes can arise due to systematic ef- 
fects. Recent detections of the cosmic shear have mea- 
sured different l evels of B modes (for a discussion see e.g. 
ilR,efregierll200l 'l. whose dominant source appears to be 
errors in the PSF corr ection. Fortunately, recent improve- 
ments in t his matter jHoekstrall207)4l lj arvis & Jam Zl 
l.larvis et alJl20p3) result in a very low level of B modes. 
Another potential systematics is th e possible intrinsic 
alignment of neighboring galaxies iBrownetall [2002; 
ICatelan et ail 1200 ll ICrittenden et aPboOli 12002^ which 
can be ameliorated by using galaxies in different redshift 
bins l|Hevmans fc Heaven J2003t iKing fc SchneiderlEool . 
or in the specific case of shear three-point fu nctions, by 
takin g advantage of its geometrical properties l|Schneiderl 
1200.1) . 



2.2. The shear two-point function 

The naive computation of the two point function of the 
shear field, ( n fijj), vanishes by symmetry reasons. This is 
because the pseudo- vector 7 transforms under rotation as 
a spin-2 object, which means that, 



7 = 7i + 172, R<t> (7) = (7i + 172) e 



i2,b 



(4) 



where represents a rotation by angle <j). It is possible, 
however, to build combinations of 7^7^ that do not aver- 
age to zero. For example, the combination {7 aIb) \AB\=e 
is invariant under rotations, and thus non-zero upon av- 
eraging. More generally, the two quantities 



{lAl*B + 1*a1b) 
{lAlB + l*Al* B ) 



AB\=8 



AB\=8 



(5) 
(6) 



are non-zero. 

This property is easy to understand in terms of a tan- 
gential and cross decomposition of the two-point shear 
correlation function. If we call <p the angle of the AB vec- 
tor relative to the x axis, we can define, for our two point 
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a tangential (+) and a cross (x) shear by 

7+ + i7x = (7i+i72) e~ i2v , (7) 

where the cross shear is oriented toward the direction 
n/A+ip, to account for the transformation properties of the 
shear. In this language the geometrical properties of the 
component of the shear are easier to understand. The tan- 
gential component 7+ is of positive parity, whereas 7 X of 
negative parity. We call parity transformation any reflex- 
ion, for example, the reflexion on the x axis that transform 
y into —y. The reflexion about the center of the segment 
AB that exchange A and B is another example that we 
will use in the following. Under parity transformation, a 
7+ will transform into itself, whereas a 7 X will acquire a 
minus sign. 

For the two-point function, a reflexion about the cen- 
ter of AB is equal to a rotation by an angle 7r due to the 
transformation properties of the shear pseudo vector un- 
der rotations. Therefore, for a parity negative two points 
function _D_ , one would have, for a parity transformation 
P 

P(D-) = -D- 

= RAD-) 
= D- 

= 0. (8) 

Only the positive parity combinations of 7+, x will be non- 
zero. In other word, any non zero two-point function can 
be decomposed as a linear combinations of 

£++(#) = (7+7+) 6 (9) 
£xx(0) = (7x7x) e . (10) 

In particular the correlation functions £± defined above 
read 

£ ± (0)=£++±£xx. (II) 

2.3. The shear three-point functions 

Here we extend the approach in the previous section to the 
three-point functions. In order to do this, we first need to 
define how we describe triplets of points. For any triplet, 
we call the three vertices 6i, with i taken modulo 3, so 
that the triangles 606162 is always oriented. We define 
the sides ii and the oriented angles at each vertex such 
that 

f-i = Oi+l — 9i+2 (12) 

^ = ■ (13) 

Figure ^ summarize these definitions. 

Because of translation invariance, the three-point func- 
tions depend only on the relative positions of the vertices, 
so we shall use the three lengths ipi to specify our oriented 
triangles. 

To define the shear three-point functions we 
decompose the shear into cross and tangential 




Figure 1. Definition and convention of triangle variables. 



components in some basis relative to the triangle 
llZaldarriaea fc Scoccimarroll2003t ISchneider k. Lombardil 
l200.1l . This is simpler than trying to find all possible 
combinations of the three pseudo-vectors 7 that have 
non-zero three-point functions, but of course is equiva- 
lent. Choosing a basis was "natural" in the case of the 
two-point function (given by the line joining the two 
points), but less obvious for a triangle. 

Since we are interested in the three-point functions, 
the knowledge of the position of the points is irrelevant; 
only their relative positions matters. Therefore, we will 
drop the vector notation for the vertices 6i. Thus, any 
configuration is completely described by either the three 
points relative positions, the three angles or the three side 
length. Of course, since we have decided to only deal with 
oriented triangles, the length, and not the side vectors are 
enough. 

A simple solution is to pick some "special" points of 
the triangle and to project the shear along the vector link- 
ing this point and the vertex we are interested in. Note 
that it is not, in fact, necessary to choose such a special 
point. Any set of three orientations defined by invariant 
properties of the configuration will do. In the following, 
we call this choice of orientation a projection convention. 
Any change of projection convention can be described by 
a rotation of each of the projection directions at each ver- 
tex of the triangle. Of course, each of these rotations can 
in principle depend on the shape of the triangle. We will 
call such change a rotation of the projection convention of 
angle {(1,(2,(3)- 

Once the projection convention is chosen, the shear 
field at the three vertices of the triangle can be de- 
composed in a tangential and a cross component, 7 = 
(7+,7x). As in the two point function case, these two 
components have, respectively, positive and negative par- 
ity. This decomposition leads to eight different three-point 
functions 7p„ p , /1, v and p being + or x . Half of them, 
7+++ and 7+ xx plus permutations are parity positive, 
whereas 7 XXX and 7x++ plus permutations are of neg- 
ative parity. Since the three-point functions also depend 
upon the projection convention, we will sometimes denote 
them by 7^S 1 p ( ' 2 ''° 3 \ (Ci, Ca> C3) being a rotation to the ref- 
erence frame. 

We have seen that for the two-point function, parity 
properties reduce the number of independent functions 
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from four to only two. One can wonder whether the parity 
properties will also reduce the number of non zero three- 
point functions. The parity negative two-point functions 
vanish because after a parity transformation, a simple ro- 
tation will bring back the points to the initial position [see 
Eq. ©]. This is not the case for three-point functions, 
except in the special cases of equilateral triangles, and 
some isoceles configuratio ns, where a rotation can mimi c 
a parity transformation l|Schneider fc Lombardil l2003h . 
This is why in the general case parity-negative three- 
point functions are non -zero l|Schneider fc Lombardil2 003: 
iTakada fc Jain1l20034 . Formally a parity transformation 
of a parity negative three-point function for a configura- 
tion 9 1 6*2 6*3 gives 



with z = 1,2, 3, which can be written as, 



P 



I ph>p 



7<£/ 2 ' C3) (4,£ 2l 4) 



= ~(Cl,C3,C2) 

/ p,pv 



(4,4, 



(14) 



where we have assumed that the transformation exchange 
the last two vertices. Thus, only if £ 3 — £ 2 , C2(^i,4,4) = 
(3(^1^2, 4) and p = parity forces the three-point func- 
tion to be zero. Note that for simple projection conven- 
tions, the property Ca(4, 4, 4) = C3 (^1,^2, 4) will usually 
be verified. 

A free parameter that we overlooked in the case of 
two-point function is the choice of projection convention. 
Indeed, we picked the easiest projection where the angles 
from the reference frame were the same for the two points. 
If we now allow these two angles to be different, we break 
the property that the (7+7x) are zero, since we cannot 
equate parity transformations and tt rotations, and we are 
then in similar situation to that of three-point functions. 
The choice of projection that seemed most natural for the 
two-point function has this extra property that its parity 
negative component is zero. 

Is it possible to find an analogous "natural" projection 
convention that reduces the number of non-zero three- 
point fun ctions? This problem has b een partly investi- 
gated by l|Schneider fe Lombardill2003T) . Here we shall re- 
produce their results and extend them to answer the ques- 
tion of the preferred projection convention. 

In the two-point function cases, as we discussed above, 
a rotation from the natural projection convention will in- 
duce non-zero parity negative part. One can show that the 
£± will transform as 



£(Cl,C2) = e -i(d±f2) 



(15) 



In other words, the £± only acquire a phase, their am- 
plitude is invariant. One can build equivalent "rotational 
inva riant" the three-point functi ons as described in details 
in l|Schneider fc Lomb ardi 2003), where it is showed that 
any change of the projections axis of the three points will 
leave the following four complex quantities invariants up 
to a phase 



To = (717273) 



Ti = 




(16) 



r„ 



A rv® — 4 F? 'V® 



(17) 



where //, v — 0,1,2,3 and the matrices A and B have all 



fl IS 



elements equal to ±1, B 
2(5 M o)- In addition, 7® and 7® denote, respectively, the 
positive- and negative-parity three-point functions, 



(18) 



The transformation of the Ti under the rotation of the 
projection directions are given by 

r = e^+^+^ro 
r', = c 2i( ~ Ci+C2+C3) ri 



/ 7+++ ^ 




^ 7xxx\ 


7+xx 


, 7® = 


7x++ 


7x + x 


7+x + 


\7xx+ / 




\7++x J 



T' 
1 



= e 2 ^ 1 



-C2+C 3 )p 



j? = c 2i(Ci+C2-C 3 ) r3 _ 



(19) 
(20) 
(21) 
(22) 



This property has some interesting consequences. Indeed, 
Eqs. I|19l22fl show that projection convention rotations 
will not change the modulus of the IV At most, if a pre- 
ferred projection exists (i.e. if the eight three-point func- 
tions are redundant), they can be reduced to four quan- 
tities. It might be possible that for a given triangle shape 
one can find a preferred convention projection for which 
all parity-negative three-point functions will be zero (one 
of course can always find a rotation that yield three van- 
ishing three-point functions). Let us note (Ci,C2,Cs) the 
rotation that will transform a given projection convention 
into the projection we are looking for. Eqs I|19l22fl define 
a system of four equations that our Q have to verify. This 
is possible if and only if $i, the phases of the Ti, verify 
the relation 



$0 = $1 + $2 + $3 + kir ke2 



(23) 



For isosceles triangles, this condition can be shown to 
be true. Indeed, the parity properties imply that for any 
isoceles configuration in 9\ , the following relations hold 



r 2 

Ft, 



1 3 

r; 
rs. 



(24) 
(25) 
(26) 



that is, $0 and $1 are zero (modulo 71"), whereas $2 = 
— $1. The rotations that change a projection convention 
to one with the zero parity negative three-point functions 
are of the form (0, — <f>/4, $/4), where $ = <f> 2 = — $3 is 
the phase in the first projection convention. Note that this 
result is of little use, since it just shows in a different lan- 
guage that only four quantities are needed to describe the 
isosceles configurations, which can be shown only from 
parity considerations. Moreover, the phase <!> has to be 
evaluated in order to find the preferred convention pro- 
jection. 

Neither parity properties, nor geometrical considera- 
tions allow us to prove or disprove the conditions in Eq. 
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l|23|) in the general case. The answer to our question will 
have to come from the computation of the three-point 
functions that we will exhibit in the next section. 

3. Analytical predictions for the three point 
functions 

We now proceed to calculate the eight shear three-point 
functions. This can be done in terms of the eight 7 M „ P or 
four I\ which are completely equivalent. In fact we will 
present results in both representations, although most of 
the time we deal with IV 

To compute the three point functions of the shear, we 
will use a simple model that captures most of the features 
of a more detailed calculation. The reason for studying 
such a simplified model will become clear in section0]when 
we discuss how to "optimally" combine the eight three- 
point functions into a single three-point function. 

3.1. A single halo 

Weak lensing surveys provide their best constraints at 
scales small enough (one to ten arc minutes) that are 
well into the non-linear regime. For example, for measure- 
ments on background galaxies at redshift around unity, 
the weak lensing efficiency is maximum at about z — 
0.43 for an f2 OT = 0.3, fl\ — 0.7 cosmology, for which 
one arc minute corresponds to a distance of 0.3 Mpc/h. 
For such scales, contributions to lensing statistics come 
mainly from light deflection by single dark matter ha- 
los, and it is a good approximation to compute the light 
deflection ignoring c o upling between different de flections 
I Coomv fc Hi J Eofll IVan Waerheke et alJ Eool. Tn the 
langu age of the so-called halo model l|Coorav fc Shethl 
2002), statistics are dominated by the "one-halo" contri- 
bution, and this has been verified for the shear three-point 
functions by co mparison with measure ments in numeri- 
cal simula tions llTaka da fc Jainll2003 aloh . R emarkably, as 
shown in 1 Zaldarriaga fc Scoccimarroll2003h . for scales of 
about three arc minutes the full dependence of the shear 
three-point functions on the triangle shape for the halo 
model agree very well with a calculation based on a sin- 
gular isothermal sphere up to an overall amplitude. Given 
these results, we will only slightly go beyond the singular 
isothermal model. We shall assume that the shear three- 
point functions can be calculated by the contribution from 
one spherical halo located at the maximum of the lensing 
efficiency window. The halo profile will be taken as a gen- 
eral power-law 



p(r) oc 



1 



(27) 



where n = 2 corresponds to an isothermal sphere. In 
general, simulations suggest dark matter halos have n ~ 
1 near the halo ce nter and n ^3 at large distances 
l|Navarro et alJll997|) . 



1 



p oc 



r(r + r) 2 ' 



(28) 



where the effective profile index is n c ff — 2 when r — 
r , typically a tenth of the halo virial radius. Assuming a 
fixed n (which can be determined by the angular scale 
probed and the halo mass of the dominant contribution 
to the three-point function) we loose cosmological infor- 
mation, but on the other hand the calculation can be done 
analytically and also the scaling properties of Eq. l(2"?|) al- 
lows us to ignore the halo profile normalization and the 
effects of projection which enter as an overall constant. 

Computing lensing by a spherical halo is very simple. 
From Eqs. © and (|2T|) it follows that the shear pat- 
tern behaves as 



7 



cos(26>) 
sin(26>) 



7(7), 



with 



7 (r) 



(29) 



(30) 



The contribution of such a halo to the three-point func- 
tions is therefore 

l^ P (hJ2J 3 ) = (31) 
J d 2 u P M • 7(00 P v ■ 7 (0 2 ) P p ■ 7(63), 

where the symbols P ^ stand for the projectors in the + 
and x directions. The three points are taken at the posi- 
tions Oi, such that ii = — #i+2, and we use all our 
indices i modulo 3. Note that we integrate out the position 
of the center of the halo, denoted by u. 

Equations (|30|> and (|31|l imply that in our simple 
model, the following scaling relationship holds, 



7^p(A4,A4,A4) = A b - J " 7 ^ p (4^2^3) 



(32) 



Since all choices of projections are equally valid, we can 
choose here the projection operator to simplify the calcu- 
lation as much as possible. Here we will take the projection 
direction at each vertex to be that of the opposite side. 
Note that this is equivalent, up to a sign, to projecting 
along the li nes joining each vertex of th e triangle to its 
orthocenter l)Schneider fc Lombardill2003h . For simplicity, 
we shall refer to our projection as "orthocenter". With 
this choice, the projected shear reads 



trig„(20i) 

P, ■ 7(0;) = - M 
16, - u\ 



(33) 



where trig + = cos and trig x = sin. The angle 4>i is the 
angle defined by the line joining the vertex and the center 
of the halo and the opposite side of the vertex 







(u - 0j) ■ tj 

\u-e,\L 



(34) 



Given our choice of projection, the trigonometric functions 
can be written in terms of the configuration (^1,^2^3)5 
and the relevant integr als computed in terms of hyperge- 
ometric functions (see ijZaldarriaga fc Scoccimarroll2003|) 
for description of the integration procedure) . We refer the 
reader to Appendix A for details. 
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3.2. Results 

Figures|3l21and0|present our basic results for the isother- 
mal profile case (n = 2), for fixed ratios of sides (£ 2 /£3 = 
1, 2, 3), as a function of angle ipi as defined in Fig.Q] Due 
to the scaling in Eq. Ij32(l a choice of sides ratio completely 
describes other triangles with different overall scale up to 
a normalization constant. Figure [21 shows 7pi, p , whereas 
Figs. |21 and 2| show T^. A comparison between different 
values of the power-law index n is presented in Figs. for 
5R(r ) and S(r 2 ) when £ 2 = £ 3 . 

Figures [21 and show results in two different projec- 
tion conventions. To go from the orthocenter to the center 
of mass projection convention it is necessary to compute 
the angles r\i between the line joining the center of mass 
with the i th vertex and the (see Fig. and use the 
relations given by Eqs. I|19I22|I . Comparing both projec- 
tions is useful to disentangle geometrical properties from 
projection-dependent behavior. Comparing top and bot- 
tom panels in Figs. [21 and [21 we see that, qualitatively, 
the orthocenter projection leads to "wigglier" correlation 
functions. 

Parity related features in the three-point functions for 
isoceles triangles are evident: 7x + x — 7xx + > 7++x — 
-7+x+ and 7 XX x = 7x++ = in fig[2]and T 2 = T% and 
3(ro) = S(ri) = in fig [21 Furthermore, for equilateral 
triangles 7 +xx = 7x + x = 7x++ and 7 X++ = 7+x+ = 
7++x = 0. Other features that these figures show appear 
more difficult to predict; for example, the local extrema of 
7+++ for configurations close to equilateral triangles are 
not exactly located at 27r/3 and depend on the projection 
convention (after all, the average over the position of the 
halo does depend on projection convention). In addition, 
points where some are equal to each other change 
location (and can disappear or appear) as projection con- 
vention is changed. 

Figure ^ shows what happens as we consider trian- 
gles other than isosceles. Note that now all the parity- 
negative three-point functions are non-zero. The fact that 
Ti ~ — T 2 (and accordingly that r ~ T 3 ) can be under- 
stood as follows. Indeed, when cos(7r — Vi) = 1/4 (1/6 
for bottom panel) the configuration is isoceles in #3, thus 
ensuring that Ti = T 2 by parity. Around this angle, it fol- 
lows from parity that T 2 (ipi + e, £ 2 , £3) ~ ^1(^1 — e, £ 2l £3) 
for e <C £2/ £3- In addition, as £2/^3 increases, the product 
7172 will start to dominate the three-point function, and 
thus parity properties become essentially those of the two- 
point function, insuring that Y\ ~ — T 2 . These arguments 
explain why as ^2/^3 increases (compare top and bottom 
panels in Fig. 0} Ti gets closer to — T 2 and T to T 3 . 

Figure shows how the correlation function 5R(Fo) and 
3(r 2 ) depend on the slope of the profile n. The zero cross- 
ing of 3(r 2 ) at 27r/3 can be explained by parity; however, 
the other zero crossings (and the number of them) depend 
on the profile slope n. For 3?(Fo) there is a zero crossing 
at ~ 50°, which appears robust to changes in n, but 
we found no simple explanation for this. 
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Figure 2. The top (bottom) panel shows the positive 
ilpup) an d negative (7® i/p ) parity three-point functions for 
the center of mass (orthocenter) projection convention. 
We have assumed an isothermal sphere profile (n = 2) 
and a fixed ratio {£ 2 jlz = 1) and plotted as a function 
of the angle tpx (see Fig. for a definition). Line styles 
for parity positive three-point functions are 7+++ (solid), 
7+ x x (long-dashed), 7 x + x (dotted) and 7xx+ (dashed) 
whereas for negative parity 7 X x x (solid) ,7 X ++ (long- 
dashed) ,7+ x + (dotted) and 7++ x (dashed). 



We now consider again the question raised at the end 
of section 12.31 regarding the existence of a preferred pro- 
jection, in the context of our simple model. The condition 
in Eq. (|2^|l can be rewritten as 

T = r (r 1 r 2 r 3 )*eM. (35) 

Figure shows 5?(~l) and 3(1) for triangle with £ 2 /£3 = 2. 
We see from the bottom panel that the imaginary part 
does not vanish, thus it is not possible to find a preferred 
projection where all negative-parity three-point functions 
vanish [but of course we can make three of them vanish 
by using Eqs. (|I9I22() ]. Note however, that the relation in 
Eq. (|35|l is close to being satisfied, at least approximately, 
for almost all ipi- 
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Figure 3. Same as Fig.FJlbut for the real (positive-parity) 
and imaginary (negative parity) parts of IV Line styles 
are as follows: To (solid), Ti (long-dashed), IV (dotted), 
r 3 (dashed). 



4. The shear three-point functions in numerical 
simulations 

We now describe the results of measuring the shear three- 
point function in N-Body simulation. For details about 
the simulations and the procedure we followed to make 
the measurements see Appendix iBl Basically, we created 
shear maps with three different resolutions that cover a 
large range of scales (roughly from 4 arc seconds up to 
1, 10 and 40 arc minutes). In practice, we measured the 
7^„p in the center of mass projection convention, and then 
transformed to the IV 

4.1. Results 

We start by comparing the scaling of the three-point func- 
tions to see, using Eq. (|32l) . whether the effective value 
of the profile index is reasonable compared to expecta- 
tions based on dark matter halo profiles such as Eq. I|28|) . 
Figure shows results from the medium resolution mea- 
surements, where we have scaled 9?(ri) assuming an n = 2 
profile. More precisely, given an isosceles triangle where 
two sides are equal to £, we have fitted for Fl and a{£) so 
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Figure 4. Same as Fig. but only in the center of mass 
projection, for sides ratio £2/^3 = 2 (top panel) and 
£2/^3 = 3 (bottom panel). 

that 

R(TifoM,*)) ~ fiWa + c^O, (36) 

where a(£) incorporates the dependence on the angle be- 
tween the two equal sides. In practice we found that vari- 
ations of a with I are small, thus they have been ne- 
glected. To avoid artificial deviations from scaling due to 
our binning (see discussion below), we have kept our tri- 
angles far enough from being collapsed, i.e. we restrict 
tt/9 < Vi < 8tt/9. The results in Fig. □ show that the 
scaling between 1 and 3 arc-min is consistent with that 
of an isothermal sphere, whereas for larger angular scales 
tT\(t) ~ £ _1 and thus the effective profile index increases 
toward n = 7/3 ~ 2.3 [see Eq. Q32[l]. This is consistent 
with the hypothesis that we are only sensitive to the in- 
dex of halos that contribute the most at the scale we are 
probing. Indeed, at the maximum of the lensing window 
function, z ~ 0.4, masses in the range 10 13 M Q to 10 15 M Q 
have ro [where n e g — 2, see Eq. (|28p] of the order of a few 
arc-minutes. 

Given the results in Fig. we can now compare the 
n = 2 model to our measurements in simulations at arc 
minute scales to check whether they agree (up to an ar- 
bitrary constant that our model does not predict). Figure 
13 shows the comparison between the n = 2 model (with 
amplitude fixed by maximizing agreement with T2) and 
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Figure 5. The real part of Tq (top) and the imaginary 
part of T2 (bottom) for different values of the profile index 
of the halo n and £2 = £3- Line styles are n = 1.2 (long- 
dashed), n = 2 (solid), n = 2.1 (dotted), n = 2.5 (dot- 
dashed) 




-0.01 1 1 1 1 1 1 1 

30 60 90 120 150 V>1 180 

Figure 6. The real (top) and imaginary part of 1 for the 
£2 1 '£3 — 2 case as a function of tpi . The non-zero imaginary 
part shows that there is no preferred projection convention 
where all negative-parity three-point functions vanish. 
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simulations for triangles with £2 — £3 — 2.25 arc-minutes. 
We see that by adjusting a single amplitude, all other I\ 
show good agreement as well, giving support to our simple 
model. This result is not surprising given that an isother- 
mal profile was shown to agree with a calculation base d 
on the halo model in ijZaldarriaga fc Scoccimarroll2003|h 
and the latter was found to be in good agreement with 
measurements in numerical simulations in l)Takada fc Jainl 



It is apparent from Fig. that as the triangle becomes 
collapsed (ipi — ► 0°, 180°) the agreement with the model is 
not as good, particularly in the parity positive case. This 
can be understood as a result of the effect of binning in 
the simulation measurements, where many configurations 
contribute to each bin. In the case we are considering here, 
it amounts to an error of order 10% on the length of £2 
and £3 and about 2% on the angle ipi . On can estimate the 
effect of binning by computing in our model I\ defined as 
instead 



p. = 1^,(1+6)1,1)+ 1^,1,(1 + 6)1) 



Figure 7. Scaling test for £q — 1.125'. The points with 
error bars represent S = £Ti(£) / (£o^i{£o)) ■ The solid line 
is S = 1 (n = 2), whereas the dotted line denotes to which 
according to Eq. I|32l) correspond to n ~« 2.3. 



with e = A£/£ ~ 0.1. This prediction is compared to the 
same measurements in Fig. ED The results show that the 
behavior near collapsed triangles, such as the sharp diver- 
gence of Sft(ri), can be explained by the discreteness error 
due to binning. Generically, our model shows that the ef- 
fect of a small error in the determination of the length of 
the triangle sides can translate into significant corrections 
to the three-point functions for nearly collapsed triangles. 
Another noticeable discrepancy between our model and 
the numerical simulations are the slightly displaced zeros 
of Q(r 2 ) and 5(r 3 ) and their amplitudes for rpi < 60°. 
The former can be improved by slightly changing the pro- 
file index, as shown in Fig. the first zero in these three- 
point functions is very sensitive to small variations in the 
profile index. The difference in amplitude can also be seen 
in similar plots in l|Takada fc Jainll20034 when using the 
full halo model; this suggests perhaps that such deviations 
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Figure 8. Predictions from our n = 2 model compared 
against measurements in simulations for isosceles triangles 
with £ = 2.25 arc minutes. The top (bottom) panel shows 
the real (imaginary) part of T±, and the model arbitrary 
amplitude has been adjusted to match T2- 



could be due to the assumption of spherical halos. It is 
worth exploring this issues further as they may provide a 
novel way of testing dark matter halo profiles and shapes. 
Some steps in this dir ection have been already taken in 
l|Takada k Jainll2003cj) . where the two- and three-point 
function of the convergence field were used to constrain 
parameters of halos such as concentration and inner pro- 
file slope. 

Finally, we present results from simulations for isosce- 
les and non-isosceles configurations, from our high reso- 
lution measurements, at £2 = 0.42' (Fig. I10[l and middle 
resolution one at £2 = 2.8' (Fig. lll|) . These should be com- 
pared with the top panel in Fig. [3] for the isosceles case, 
and Fig. 0] for non- isosceles triangles. We can see that for 
£2 = 2.8' the agreement with Figs. I314I as expected from 
the scaling test that suggests that the effective index is 
close to n = 2. However, for £2 = 0.42' we see that there are 
significant differences, the three-point functions changed 
as expected from the results of our analytical model in 
Fig. when the profile index becomes smaller than n = 2. 

We have also made measurements in lower resolution 
sets (not presented here) that allows us to probe larger 
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Figure 9. Same as previous Figure, but modifying the 
model prediction to mimic the effect of binning in the 
simulations, approximately, as in Eq. 1371) . 



scales. Again we see consistent results with those expected 
from Fig. [S] for indices n > 2, in particular regarding the 
dependence of number of zeros as the scale is changed. 
However, one cannot extend this study to significantly 
larger scales, as contributions from more than a single halo 
become important and our simple model breaks down. 

5. Estimators of cosmic shear three-point 
functions 

We now turn to the problem of building a simpler esti- 
mator of the shear three-point function. One obvious so- 
lution is to combine the m to reconstruct the w eak lensing 
convergence bispectrum ijSchneider et al.ll2005|) . However, 
equations involved into the computation of the bispectrum 
from the shear three-point functions are very difficult to 
compute numerically. They correspond to the inversion of 
a non-local equation; a task that usually cannot be ful- 
filled from real data without some regularization of the 
problem. 

The so-called M ap statistic filters provide such regular- 
ization. They have been used successfully to reconstruct 
the filtered convergence two-point functio n from measure- 
ment of the shear two-point functions l|Schneider et alJ 
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Figure 10. I\ measured for £2 = 0.42'. In each panel top 
(bottom) part shows positive (negative) parity three-point 
functions. Panels correspond to isoceles triangles (top), 
£2 — 2^3 (middle) and £2 = $£3 (bottom). Symbols are as 
follows Tq (plus),!?! (cross), T2 (star) and T3 (box). 



2002]). Another great advantage of these family of com- 
pensated filters is that, to some extent, they can be cus- 
tom designed to have a finite real space support (or at 
least exponentially small wings), allowing for relatively 
quick data analysis. A problem of those method, however, 
is that one can expect a degraded signal to noise from the 
initial data, due to the use of compensated filters that im- 
pose cancellation of part of the signal. This have not been 
an important issue for two-point functions. 



Figure 11. Same as previous Figure for £2 = 2.8'. 

The M ap appr oach has be e n als o applied to three 
point functions lljarvis et al.l 12004 IPen et all [2003; 
ISchneider et al.ll2005|) . It has been showed that one can 
exhibit a summation procedure for measured data that 
results into an estimation of the M ap filtered bispec- 
trum of the convergence field. Measurements on real data 
have been performed. The quality of data being low and 
the degradation of the signal-to- noise ratio inherent to 
M nri statistics results in very large error-bars l|jarvis et alJ 
l2004UPen et"alll2003l) . 

Owing to its simple relation with the convergence bis- 
pectrum and thus to the matter distribution bispectrum, 
the M ap approach is certainly a very appealing estima- 
tor of the weak lensing three-point function. However, as 
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it seems to require a high quality dataset, the need for 
a simple and robust estimator of the weak-lensing three- 
point function remains. The shear three-point functions Tj 
can provide such tool, but they still have a complicated 
dependency on the configuration. Our goal in this section, 
is to use our analytical prescription in order to propose 
a way to build an estimator with simpler properties, yet 
avoiding as much as possible cancellations to preserve the 
signal-to-noise ratio. 

5.1. Previous estimators 

Bernardeau, van Waerbeke and Mellier l|Bernardeau et all 
l2003ft (hereafter BvWM) proposed a simple estimator that 
exhibits some of the properties discussed above. They 
studied the pseudo vector field T 

= <7i (7 2 -73))= (38) 

where all pseudo vectors 7 are projected along the direc- 
tion given by t\. In such case, symmetry does not im- 
ply that T vanishes. BvWM were able to evaluate T for 
special configurations of the three points. They also com- 
puted it in the framework of the hierarchical ansatz and 
measured it in N-body simulations. From their results, it 
seemed that the parity-negative part of T was negligible 
and that the parity-positive part was not changing sign 
in a small region around 82 83. They empirically fitted the 
shape of this region by an ellipse of focal points 82, 83. 
In light of this, they proposed to use the averaging on 
the ellipse of the parity positive component of T as their 
estimator of the shear three point function. 

The pseudo-vector F can be re-expressed as a combi- 
nation of the IV Let us place ourselves in the projection 
convention defined by the direction of 8283 (this choice 
of projection is peculiar in the sense that it artificially 
breaks the symmetry properties by distinguishing one of 
the points). With this choice, the scalar product (j 2 • I3) 
reads 

(72 ' 73) = 7(2) + 7(3)+ + 7(2)x7(3)x- (39) 

Multiplying the previous Eq. by 7! leads to the expression 
of T 

T= (7+ ++ +7 + xx \ (40) 
V7x++ +7xxx / 

which can be written in terms of the I\ 

■F = i(r 2 + r 3 ). (4i) 

Remember that the last equation holds only for the 8283 
projection convention; Eqs. i|19l22[l can be used to rotate 
the expression into some other projection convention. 

FigureElP resen t s the T pseudo- vector in the top right 
quadran of the plane, computed from our analytic model, 
with index n = 2. The base of triangle (82 6*3) spans the 
range [—1/2, 1/2] on the x axis. Fig. IT2l p artly reproduce 
the result from l|Bernardeau et alJ l2003j) . In particular, 
even if the parity-negative part of the pseudo-vector T is 
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Figure 12. Predictions for the T estimator. Left panel 
presents the real/parity positive (top) and imagi- 
nary/parity negative (bottom) parts. Right panel shows 
the resulting pattern. The points 82 and 83 are fixed to 
(0,-1/2) and (0,1/2). The coordinate gives the position 
of 8\. We only shows one fourth of the plane. The rest of 
it can be determined from parity properties. T x goes go 
to zero on the x and y axis. 



smaller than the parity positive component, it is not negli- 
gible. We also explore the region where the parity positive 
component keeps the same sign. Contour plots on the left 
panels of Fig. ^Jshows that its shape is richer than that of 
an ellipse. Changing the profile index n modifies slightly 
the properties of T . The overall shape shown Fig. ^| is 
conserved. Modifications are concentrated around the x 
axis, when the configuration is nearly flat; in particular, 
the number of constant sign region along yg 1 = changes 
as can be expected from Fig. [5] Finally, the amplitude of 
the parity-negative component decreases slightly when the 
value of the profile index increases. This explains the ap- 
parent discrepancy between BvWM results and ours. They 
are averaging their measurement in numerical simulations 
on scales bigger that the one probing the region where the 
halo profile is well described by n = 2. Moreover, the bin- 
ning procedure can average out the parity negative part. 

Although ignoring the parity negative part and sum- 
ming the parity positive part over an ellipse is a good 
starting point, the BvWM estimator can be improved by 
using the prediction of the expected shear pattern as we 
now discuss. 

5.2. A new estimator: projecting measurements onto 
shear pattern templates 

An approach to avoid signal to noise cancellation due to 
positive and negative contributions when combining differ- 
ent three-point functions is to "project" the measurements 
of the r, directly onto the expected result (template) from 
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analytic predictions and build an estimator such as 

p = ^f,r*, (42) 

i 

r, being the analytical predictions, the star denoting the 
complex conjugate. Note that T>q now depends on our 
choice of profile index n. One can do even better, assum- 
ing that the covariance of the Ti is known, by building a 
"minimum variance" estimator of the three-point function 

x^^r^r*} -1 ^*, (43) 

where ^r,T*) 1 denotes the inverse of the covariance ma- 
trix. T>\ is only a minimum variance estimator to the 
extent that its distribution can be approximated to be 
Gaussian. 

It has been shown that to a good approximation, the 
covariance of the three-point function of the weak lens- 
ing shear can be evaluated by re stricting the computa - 
tion to the Gaussian contribution IjTakada fc Jamll2003q) . 
However, even if we restrict ourselves to the Gaussian con- 
tribution, the covariance matrix is difficult to compute 
analytically. Indeed, as discussed above the geometrical 
properties of the shear complicate the calculation of the 
three-point functions, and the situation is even worse here, 
as we have to integrate over all possible orientations and 
positions of two identical triangles. To avoid such compli- 
cation, we evaluate the covariance matrix by measuring 
it in our numerical simulations (using 40 realizations, see 
Appendix B). Even in this case the resulting evaluation of 
the covariance matrix is quite noisy, which somewhat re- 
duces the reliability of 2?i, but nevertheless provides an 
estimation of the expected signal-to-noise improvement 
that can be expected from a better determination of the 
covariance of the IV 

The T>o and T>\ estimators depend on the configura- 
tion of the three points. We reduce this dependence by 
summing the estimators over a set of configurations. This 
sum is similar to the approach taken by BvWM, where 
they integrated over all observed configurations in a small 
ellipse. In our case there is no particular choice of summa- 
tion region, as we are guaranteed to avoid cancellations as 
long as the effective profile index is close to the one of our 
template (n = 2, our fiducial choice). To simplify the pro- 
cess of summation, we take advantage of the fact that our 
measurements are already binned by length and opening 
angle, thus we sum configurations along the opening angle 
of the triangle, ipi, for a given ratio of lengths £2/ £3- The 
overall length of the triangle should be absorbed into the 
scaling relation (|32[1 when it holds. Thus we define 

zb(W4») = / <tyi v (tpi,e,ex £ 2 /4), (44) 

and the corresponding X\ , for the estimator T)\ . When the 
effective profile index is close to n = 2, Eq. I|32|) implies 
that To(£ , £2 1 '£3) obeys the following relation 

2o(AW4) = 4 X o(W4)- (45) 



In other words, we expect the estimators to behave like 
power laws. Figure 1131 and 1141 present the measurements 
of 2o an( l X\ with the template profile index n = 2. Error 
bars correspond to the variance among different realiza- 
tions. The integral over the opening angle has been cut- 
off at tpi < 7r/9 and Vi > 8tt/9 to avoid errors aris- 
ing from binning. The measurements are compared with 
analytic predictions obtained by mimicking the measure- 
ment procedure, i.e. the integrals of V^ 1 — J2i ^i^i (resp. 

vt \ = E l j ( r i r j) _1 T*) are computed with the same 
cutoff on ipi and using the same sampling in ip\ resulting 
from the binning of our dataset. This explains the slight 
departure from a power-law in the theoretical curves in 
Fig. 1131 The normalization free parameter in the theoret- 
ical model has been set so as to maximize the agreement 
between the measured Y2 and its analytic estimation for 
isoceles configuration at 1 arc-min. The large difference 
in normalizations between X an d ^1 is due to the fact 
that we did not normalize the covariance matrix; it ac- 
counts for the inverse of the covariance matrix determi- 
nant. The estimators Xq and X\ have been measured for 
isoceles triangles, as well as for some elongated configura- 
tions £2^3 = 2, 3 and 4. Due to our data analysis strategy 
only a few elongated configurations are available (see dis- 
cussion in Appendix [EJ) . Figure ^5] shows the results for 
all our dataset, whereas Figure [21 focuses only on the 
medium scale dataset, where the effective profile index is 
expected to be close to n = 2. 

As expected, Iq is very close to its theoretical estima- 
tion between 1' and 4' where the effective profile index 
is close to n = 2. The agreement is very good for isoce- 
les triangles as well as for the elongated configurations. 
The agreement between the theoretical and measured X\ 
is not as good. Note that Figure [21 seems nevertheless to 
indicate a behavior similar to a power law, corresponding 
to n c ff — 2. Clearly, here we are sensitive to the noise in 
our estimation of the covariance matrix. Even with our 
cut in ipi , T>i still gets a significant contribution from flat- 
tened configurations, where the covariance matrix terms 
are large. The discreetness errors induced by the data 
binning are also responsible for part of the discrepancy. 
Measurements for elongated configurations {£2/ £3 > 1), 
where discreetness effects are less important by construc- 
tion, show a better agreement with the analytic estima- 
tions. For scales out of the 1 — 4 arc-min ranges, the the- 
oretical pattern is no longer valid, and the projections T>o 
and T>i decrease. We can see this behavior in Figure H~3l 

We now estimate the improvement of signal to noise in 
our estimator compared to the one proposed by BvWM. 
More precisely, we measure the pseudo-vector J 7 , as a 
function of the configuration, project it on our analytical 
model and sum it along the opening angle of the trian- 
gle, in a similar way to what we have done for Xq and X\ . 
Note that this is an enhancement of BvWM method since 
we are using both components of the pseudo vector T . 
Indeed, since we are projecting on analytical predictions, 
the parity-negative part will no longer average to zero. 
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Figure 13. 2q measured in our three datasets. Results for 
all datasets are combined into this composite figure. In the 
top panel, top to bottom set of points corresponds to mea- 
surement for £2 — r£$ for r = 1,2,3,4. The corresponding 
theoretical prediction for n = 2 are shown on top of the 
points for r = 1 (solid line), r — 2 (long dashed), r = 3 
(dashed) and r = 4 (dotted). Bottom panel shows the ra- 
tio between the measurement and theoretical prediction 
for r = 1. 

Figure 1151 presents the measurements in the simulation 
as well as the theoretical prediction and can be directly 
compared with Figure H31 

Finally, Figure QUI shows the signal to noise ratio for 
the different estimators, evaluated in our simulations. We 
do not take into account experimental noise such as the 
distribution of intrinsic ellipticity of the galaxies or the 
uneven distribution of sources. Our estimation is thus 
dominated by the cosmic variance. The shot noise term 
due to the intrinsic orientation of the galaxies will mainly 
contribute as a non correlated source of noise at small 
scales. For the simplified case of the convergence three- 
point function of equilateral triangles, it has been shown 
that thi s term contributes at scales smaller than 1 arc- 
minute ijTakada fc JainlfeoOSd) . As expected, X\ has the 
best S /N ratio, nearly twice better than Zq . This however 
degrades quickly as the effective profile index leaves the 
n e g ~ 2 region. The S/N ratio of Zq is about 30% better 
than for our improved BvWM estimator. This is not un- 
expected, Eqs. (|41|l shows that T can be computed with 
only two of the whereas Iq uses all four of them. If I\ 
are uncorrelated, this should correspond to a \[2 degrada- 
tion of the signal-to noise ratio, which is close to the 30% 
we measure here. 

Recall that our improved BvWM estimator uses the 
full T pseudo- vector. Using only the parity positive coor- 
dinate will further reduce the signal to noise by another 
factor y2 if its two components are uncorrelated. 

6. Conclusion and Discussion 

We have described a simplified analytic model that can 
be used to compute the three-point functions of the shear. 
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Figure 14. X\ measured in our medium resolution 
dataset. The solid line is a power law proportional to l/£ 2 
and shows the expected scaling. Other curves are semi an- 
alytical results (the covariance matrix has been evaluated 
from the simulations). They correspond to r = 1 (long 
dashed), r = 2 (dashed), 3 (dotted), r = 4 (dot-dashed). 
Measurements partially agree with the semi-analytical 
predictions 
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Figure 15. Same as Fig. 1131 for our improved version of 
the BvWM estimator. 



This model is inspired by halo models, and only considers 
the one-halo contribution of a spherical potential of power 
law profile. The free parameters of our models are the 
profile index n and a normalization. We have used this 
model to investigate some geometrical properties of the 
shear three-point function. In particular, we have shown 
that there is no preferred projection choice that will reduce 
the number of independent three-point functions. 

We have compared the predictions of this model with 
results from N-body simulations. The predicted and mea- 
sured three-point functions of the shear where shown to 
be in good agreement. In particular, we have shown that 
the approximation made on the profile index is sufficient 
to predict the shear three points function in a reasonable 
range of scales. The isothermal profile case, n = 2, corre- 
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Figure 16. Signal to noise for Tq (solid line), T\ (dashed) 
and our improved version of the BvWM estimator (dot- 
ted) . Only isoceles configurations have been taken into ac- 
count. Data are projected on the n = 2 model predictions. 
Only the medium resolution datasets are used. 



sponding to scales from 1' to 4'. These scales correspond 
indeed to the scales where the dominant halos at redshift 
z — 0.4 are seen with a local profile index n — 2. As ex- 
pected, the smaller scales exhibited a behavior compatible 
with an index n < 2, while larger scales where compatible 
with n > 2. Our model allows for computations with a 
modified profile index, and it can thus be used at smaller 
or larger scales. However, at scales bigger than 4', the one- 
halo dominant contribution model will break down, and 
one should take into account two and three halo terms 



(Takada & Jai 



Using the agreement between our model and synthetic 
data, we proposed an optimized measurement of the cos- 
mic shear three-point function. Our method trades the 
easier cosmological analysis allowed by M ap reconstruc- 
tions methods for a better signal to noise of the measure- 
ment by avoiding cancellations. We use the model predic- 
tions to compute "optimal" weighted sums of the eight 
three-point functions. Contrarily to the M ap statistic, this 
kind of estimator is local and is not affected by the shape 
of the survey. These methods can be seen as a refined 
versio n of the one proposed by BvWM l|Bernardeau et al] 
1200.1) . 

We computed two such estimators, the first a sim- 
ple projection on our theoretical model, the other tak- 
ing into account an estimation of the covariance matrix 
of the shear three-point functions. We compared them 
with an improved version of the estimator implemented 
by BvWM. Our estimators perform much better than the 
improved BvWM. The minimal variance estimator lets us 
expect more than a factor of two gain in the signal to noise 
in the best case. 

With future space-based experiments, the quality of 
cosmic shear data will greatly increase and the loss of 
signal to noise inherent to compensated filters will be a 
not too high a price to pay for accessing to the filtered 



three-point function of the convergence. In the meantime, 
we believe that methods using projections of data onto 
theoretical predictions, as the one we proposed here, will 
be an interesting alternative. Improvement to what we 
proposed here will be doubtless needed. In our analysis, 
we first measured the three-point functions and then pro- 
jected them on analytic templates, increasing the errors 
due to discreteness of the binning. We saw that this error 
has a significant impact on elongated isoceles configura- 
tions. Projecting the data as they are measured can solve 
this problem. Using a simple projection already improves 
the situation compared to what has been done before. The 
minimal variance estimator promises an even better abil- 
ity to detect the shear three-point functions. Results from 
this estimator are yet difficult to forecast as our evaluation 
of the covariance matrix is somewhat noisy. More numer- 
ical simulations will be needed to quantify this better. 

Improving the choice of effective halo profile index will 
also result in a net improvement, since by restricting the 
index to that of an isothermal sphere, n = 2, we were only 
able to improve the efficiency of the three point function 
estimator around 1-4 arcminutes. This is not a restric- 
tion of the model, since it can be easily extended to any 
other index profile. More generally, given any model for 
the three-point functions one can compute the "optimal" 
estimators we define here. 

Finally, we have not investigated how our method can 
be used to obtain information on the underlying cosmo- 
logical model. With our simple model, all the cosmological 
information is encoded in the free normalization param- 
eter. Further work, comparing our simple model, with a 
full halo model will be needed to investigate this point. 
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Appendix A: Projection onto the opposite side 

Here we describe the computations of the trig M (2<^) fac- 
tor for the opposite side projection described in section 3.1. 
Remember that the angles <j>i are defined by 



cos(^)i) 



\u-0i\h 



The computation of trig + (2<^) = cos(2^) is the simplest. 
Defining di = 6i — u, one gets 



COS((pi) 



Oi ■ t-i — 5 (#i+l — — $i+2 + d i+ 2) 



which, using the identity 

29i ■ £i — Qf +1 + Ot+z = cj + i — i i+ 2 
reduces to 

j2 j2 , »2 p2 

t i n a i + l ~ a i+2 T <-j+l — *i+2 

C ° SW = ^ ' 
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The last step needed to obtain cos(20i) is straightforward, us- 
ing cos(20i) = 2cos(</>i) 2 — 1, one gets the final expression of 
cos(20j) as a polynomial of the 4th degree in d i+ i and di+ 2 
and second degree in di, divided by d 2 £ 2 



trig+(20i) = K+i + rf *+2 ~ 2d 2 +1 d 2 + 

+2 (d 2 +1 + d 2 +2 ) (£ 2 +1 - e i+2 \ 



2d 2 



The parity negative terms require a little more work. One 
can use the identity sin(2cjf>i) = 2 cos(0i) sin(^). The difficulties 
lie in the computation of sin((/>i). One can use the fact that 



. ,, . {u-0i)A£i 

and re-express the unit vector z as 

£i+2 A ti 



(A.l) 



(A.2) 



£i+2ii sin(i/)i+i) 
Finally, a careful invocation of the identity 
(A A B) ■ (C A D) + (B ■ B) (C ■ D) = (A.3) 

(A A C) ■ (B A D) + (A ■ C) (B ■ D) 
allows us to re-express Eq. JA. Ill as 

Sin(0 ° = d^ +2 sin(^ +1 ) 

+ [(it - 6> l+2 ) + (0 i+2 - Oi)] ■ i i+ 2 1 2 } (A.4) 

which can be written as a function of cos(0i+2) and cos( 
and a term that only depends on the configuration 

1 2 

[COs((f>i +2 )di +2 £i +2 £i 



«2l 



dd 2 £ l+ 2 sin(?/>i+i) 
+ cos{</>i)di£i £i ■ ir+2 + £i+i ■ ti+2%) ■ (A.5) 

This last expression can be expressed only in terms of £7=1,2,3 
and dj=i,2,3 as the ratio between a 4th degree polynomial in 
the dj = i,2,3 divided by di, the coefficient being functions of the 
£7=1,2,3- After simplifications, the final result for trig x (20i) = 
sin(20i)reads 



trig x (2<^ 



{d 2 +i — d 2 +2 + t 2 +i ~ $+2) 
2d 2 £ 2 ^£ 2 +l £ 2 +2 -{£ 2 -£ 2 +1 -£ 2 +2 ) 2 
x [— c ld i £ i + (£ t + — ^i+2) 

+d 2 +2 (4 - el+i + £+2) + 4+i + 4+2 

-£ 2 {£ 2 +1 + £ 2 +2 ) + 2£ 2 +1 £ 2 +2 ] (A.6) 



Appendix B: Numerical simulations 

alll200lft code to evolve 24 



We used the GADGET (|Springel 
realizations of the large scale structure in a small section of a 
FRLW universe with parameters £l m — 0.3, £Ia = 0.7 h — 0.7. 
The initial conditions were imposed at redshift z = 50 us- 
ing se cond-order Lagrangian perturbation theory dScoccimarrd 
1998). The initial power spectru m used was obtained using the 
B&E jBond fc Efstathioull984T) fitting formula and normalized 
to as = 0.9 at 2 = in linear theory. The boxes we consid- 
ered are cubes of (100Mpc//i) 3 containing 200 3 particles. This 
choice has been made has a trade-off between speed and accu- 
racy. The simulation were run on 24 nodes of our cluster. For 
each realization, we took a snapshot of the large scale struc- 
tures every 100 Mpc//i. 



B.l. Building the light-cone 

We use these snapshots to compute the weak lensing effect at 
a redshift of unity on a square light cone. At z = 1 the sides 
of the boxes represent 2.47 degrees. We will only produce 2.4 2 
square-degrees patches of the sky. 

The line of sight to the sources, is built by tilling 17 snap- 
shots at different redshifts. We use this construction to create 
more pseudo realizations of the lensing effect that we have dif- 
ferent realizations of the density field. We follow a very similar 
method to the one exposed by White and Hu in JWhite fc Hul 
1200011 . 

While building each line of sight, we will make sure that 
each N-body realization is only used once. To increase the ran- 
domness of our lensing pseudo-realization, we translate each 
snapshot by a random vector, taking advantage of the periodic 
boundary condition of the boxes. Moreover, we trace the path 
of the photon along a random direction in the box, and not 
along the axis direction. Note however that all rotation angles 
are not allowed if one wants to avoid tracing twice the same 
structures in a given box. We take this problem into account 
when picking the ray-tracing direction. 

The projected mass density is then built as follows. For 
each of our 17 snapshots, after translation and rotation, we 
build a list of particles belonging to the light cone. This list of 
particles is flattened onto a 2D density map, using a Cloud-in- 
Cell algorithm. This 2D map is then multiplied by the efficiency 
window function of the lensing effect to provide a map of the 
convergence n. These k slices are added to produce the final 
lens effect on the source plane. The shear field 7 on each of 
these synthetic k maps is obtained by num erically solving Eq. 
© with a FFT llVan Waerbeke et alJlioOll) . 

Note that in our computation of the lens effect we neglected 
some secondary effects. For example our ray-tracing scheme 
is strictly restricted to the Born approximation. We thus as- 
sume here that the lens effect computed along the unperturbed 
path of the photon gives a good evaluati on of the effect. This 
as been tested in num erous work before llVan Waerbeke et alJ 
l200ll : |jain et al]l2000|) . Doing this, we neglect the lens lens- 
coupling which is known to produce non-zero corrections to 
the three-point function of th e convergence field. This c orrec- 
tion is expected to be small (I Van Waerbeke et alJr200ll) . and 
we will neglect it here as our goal is mainly to evaluate the 
validity of our simplified halo model. 

The number of slices of our light-cone is more of a concern. 
A naive evaluation of the impact of this choice can be made 
by comparing a step summation of the filtered lensing power 
spectrum with its full integration. For an Einstein-deSitter uni- 
verse, assuming a power law matter density fluctuation, we 
thus have to compare step summation and integration of the 
function 



R(t) =x 4 - n (x~lj 2 



(B.l) 



For the spectral index n ~ — |, the difference between the 
summation and integration falls below one percent as soon 
as the number of sl ices is bigger than 5. Vale and White 
llVale fc White! 120041) recently performed a less naive eval- 
uation of the same effect using N-body simulations. They 
showed that with simulations comparable to ours (300Mpc 
and 512 3 particles) one can expect about 5% discrepancy on 
the power spectrum between a computation with a slice every 
As = 125Mpc (we are slicing every lOOMpc ) and one with 
As = 25Mpc. 
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Figure B.l. Power spectrum of k measured on the 40 
lensing pseudo-realizati on. Lines are the predic ted power 
spectrum using the PD jPeacock fc Doddslll996l) (dashed) 
or halofit IjSmith et al]l2003fl (solid) prescription for the 
small scales. Points are the measurements in the simu- 
lations, without smoothing (crosses), and with smooth- 
ing (xs) to remove the Poisson noise (see text). Bottom 
panel shows the ratio between the smoothed measured 
power spectrum and the two predictions using PD (xs) 
and halofit (crosses). 



B.2. Resolution 

During the Cloud-in-Cell remapping of the particle list, the 
interpolation was done on a 2048 2 grid. The maximum reso- 
lution of our simulation is thus 4.21 arc-sec. However, at this 
scale we expect to probe the regi on where shot no ise starts to 
dominate [see equation (24) from lljain et alJl200fi) l . To reduce 
shot noise contributions we smooth the k maps by a 2 pixel 
wide Gaussian window. 

The average over our 40 realizations of the k power spec- 
trum is presented Fig. IB.ll and compared with the theoretical 
one obtained by two ansatz of the non-linear density power 
spectrum. The power spectrum is somewhat higher than the 
analytical predictions in the non-linear regime. A similar be- 
havior can be seen in the 3D power spectrum of our simula- 
tions. 



B.3. Measurement strategy 

The measurement of two and three point functions in real space 
requires significant computing time. These operations scale re- 
spectively as N 2 and N 3 where N = 2048 2 in our case. To 
reduce the time needed to perform the computation, we will 
only probe configurations up to a cut-off scale. This reduces the 
amount of computer time needed to perform the measurement, 
but prevents us to access to a broad range of scales. 

To reduce the computation cost yet preserving the abil- 
ity to measure the three-point functions over a large range of 
scales, we measure the two- and three-point functions at large 
scales on degraded resolution maps obtained by top-hat fil- 
tering and regridding the synthetic shear fields. For each field, 
we produced three datasets; one with the nominal resolution of 





max 


medium 


min 


A£ 


4.22" 


16.9" 


1.125' 


^max 


67.5" 


9.0' 


36' 



Table B.l. Description of the three datasets. Max is the 
best resolution dataset, min the worst one. Thorough this 
article, we will mainly use the medium resolution. 



2048 2 pixels representing of 4.22 2 arc-sec 2 , one with a four time 
degraded resolution (512 2 pixels, 16. 9 2 arc-sec 2 ) and the last 
with a sixteen time degraded resolution (128 2 pixels, 1.125 2 
arc-min 2 ). For the two-point functions, the computation cost 
is lower and we are able, for the two last resolution sets to 
measure it without cut-off scale; we apply a cut-off only for 
the biggest map and only consider scales smaller than a 32th 
of the map size. 

For the three-point functions we only explore a small region 
of the three points configuration space. Table |B~T1 presents the 
scales probed by each of our datasets. 

The cut-off scales have been chosen so as to have each mea- 
surement to require about ten days of computation time of one 
node of the our cluster. 

We take into account the cut-off scale we added to op- 
timize the measurement. Since the data are gridded, we can 
precompute for any given point all the positions of the couples 
of points which will create a valid configuration. By valid con- 
figuration we mean any triangle that fits the requirement of 
the cut-off scale, whose three point function measured with a 
given projection convention has a meaningful result (i.e. for the 
center of mass projection, we throw out configurations where 
the center of mass is one of the vertices of the triangle) and is 
such that a given configuration is only seen once when varying 
the initial position. This last requirement can be fulfilled by 
requiring that for any couple of points 

xi > Xi, x 2 > xi yi > yi if x\ = Xi, 
y > yi if X2 = xi, 



(B.2) 

provided that the initial point will go through the grid increas- 
ing its yi position first. 

We build a table containing for each of such configuration, 
the offset of the points positions, as well as the projector vec- 
tors on the + and x direction for each point of the triangle. 
Populating this table is an expensive task as it goes as I 4 in 
time and memory, where I is our cut-off scale. Once this initial- 
ization done, however, we just have to traverse the shear map 
and apply for each point the rules contained in the initializa- 
tion table. 



B.4. Shear two-point functions 

We measured the two point functions of the shear £± 

£±(0) = <7+7+(0)>±<7x7x(0)>, (B.3) 

where the shear pseudo-vector is projected to the + and x 
directions defined by the vector linking the two points and 
the same vector rotated by 7r/4 (see section l2~2ll . The relation 
between the shear two-point funct ions and the power spectrum 
of the convergence is well known jKaiseJll992T) 



/ 



&UP K {1)J OA (10). 



(B.4) 
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Figure B.2. The shear two-point functions in the simula- 
tions. Lines are the predicted power spectrum, taking into 
account the survey size and resolution. Solid line is £+, 
dashed, Points show the measurements in the 40 real- 
izations, for our three different datasets. The open circle 
(resp. plain squares, open squares) shows £_ and regular 
pluses (resp crosses, dashed pluses) £ + , measured in the 
max (resp medium, min) resolution datasets. 



Figure IB.2I presents the measurement of f± in our sim- 
ulations as well as analytical predictions. It can be surpris- 
ing to see that at a scale corresponding to only a fourth of 
the box scale, the two point shear correlations, and especially 
the one, differ greatly from the usual analytical prediction. 
This corresponds to the effect of the finite survey size. Even 
at a scale four times smaller than our survey size, we are al- 
ready suffering from the lack of large-scale correlations. We 
can reproduce it in our analytical prediction as follows: we are 
filtering the convergence power spectrum by a Bessel Jo or 
J4 function. Jo and J4 are maximum when their argument is 
small. This means that when the angular scale is large, the 
filter overweights the small I in P(l). This is more severe for Jo 
as its envelope decrease more quickly than the one of J4. Even 
though I P(l) quickly decrease when I goes to 0, it amounts to 
a difference that can be observed in the comparison between 
analytical computations and measurements. If we artificially 
cut the analytical power spectrum in our computation to re- 
produce the absence of scale larger than y/2 x 2.4 arc-minutes 
in our simulation, we obtain a far better agreement with our 
measurement. 

Similarly, the small scale behavior of £± is modified by 
the filtering we applied to our synthetic field by downgrading 
their resolution. The analytic predictions, once this filtering is 
included, are in good agreement with our numerical results. 

We are thus quite confident than our implementation of 
the measurement is sound and that the resolution degradation 
procedure gives meaningful results. 
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